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We calculate thermal transport in the Falicov-Kimball model on an infinite-coordination-number 
Bethe lattice. We perform numerical calculations of the thermoelectric characteristics and concen- 
trate on finding materials parameters for which the electronic thermoelectric figure-of-merit ZT is 
large, suggesting potential cooling and power generation applications. Surprisingly, the Bethe lat- 
tice has significant qualitative and quantitative differences with the previously studied hypercubic 
lattice. At low temperature it is unlikely that these systems can be employed in thermoelectric 
devices due to the low conductivities and due to a larger lattice contribution to the thermal con- 
ductivity kl, but at high temperature, the thermoelectric parameters appear more promising for 
devices due to a significant enhancement of ZT and a smaller relative contribution by the lattice 
thermal conductivity. 
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Materials that are candidates for thermoelectric 
applications 1 are interesting because they have the po- 
tential to compete with conventional compressor-based 
coolers. These materials are characterized by the ther- 
moelectric figure-of-merit ZT, which evaluates the effi- 
ciency of a bulk material as a cooling element. Com- 
mercially available semiconductor devices exhibit a ZT 
of about unity. Such values of ZT are far too low to be 
competitive with the efficiency of conventional mechani- 
cal refrigerators, which would require a ZT of 3 or 4. 

Current directions of research include bulk semi- 
conductor devices 2 ' 3 , semiconductor nanostructures 4 , 
metal/correlated semiconductor heterostructures J , and 
strongly correlated materials 6-9 . Bulk semiconduc- 
tor devices are most commonly employed, and ZTs 
of up to 1.14 have been achieved for a p-type 
(Bi2Te3)o.25(Sb2Te3)o.72(Sb 2 Se3)o.o3 alloy at room tem- 
perature and pressure 2 . Below room temperature, how- 
ever, a maximum ZT of only 0.8 has been found in bulk 3 
CsBi4Te6 at 225 K. Semiconductor nanostructures ap- 
pear more promising by increasing ZT due to a sup- 
pressed lattice contribution to the thermal conductiv- 
ity. The value of 2.4 has been observed experimentally 4 
for p-type Bi 2 Te 3 /Sb 2 Te3 superlattice structures at 300 
K. These structures also showed improved values of 
ZT (~ 1.7) at temperatures as low as 210 K. Never- 
theless, at lower temperatures all measured bulk ma- 
terials exhibit ZTs that are lower than unity. A re- 
cent theoretical proposal predicts an enhanced ZT in a 
metal/correlated semiconductor structure 5 . Estimated 
values of the low temperature ZTs are enormous; for 
instance, for a semiconductor gap of 100 meV, the best 
value of ZT would be ~ 6 at 150 K and - 100 at 40 K. Fi- 
nally, strongly correlated materials have been suggested 
as possible candidates for bulk thermoelectric materials, 
especially at low temperatures 6-9 because the correla- 
tions can strongly renormalize the Fermi temperature to 
low temperatures. Very recently, Peltier cooling has been 
achieved experimentally 9 below 10 K using crystals of 



the Kondo metal CeBg. The greatest value of cooling 
at T = 4.5 K was equal to 0.2 K, corresponding to a 
ZT of about 0.26. However, the properties of such cor- 
related materials are still not well characterized, namely, 
whether large values of ZT predicted theoretically can 
be realized experimentally or if they are the result of un- 
avoidable simplifications. 

In this paper, we discuss electronic (thermal and 
charge) transport for the Falicov-Kimball model on an 
infinite-coordination-number Bethe lattice. In a re- 
cent article 8 , we investigated transport on an infinite- 
dimensional hypercubic lattice and found possible sets of 
parameters where ZT > 1 for a wide range of tempera- 
tures. However, some of the Bethe lattice's properties are 
closer to those of real three-dimensional systems. First of 
all, the electronic band structure for the Bethe lattice has 
a finite bandwidth, as opposed to the hypercubic lattice 
where the interacting density of states (DOS) decreases 
exponentially. Another important distinction is that the 
"quasiparticle" scattering time t(u>) for the Bethe lat- 
tice is exactly zero whenever the interacting DOS is zero. 
This is also expected for a real physical system, in con- 
trast to the hypercubic lattice 8 where the scattering time 
approaches a non-zero constant as w — > ±oo and behaves 
as a power law at small oj in the "pseudo gap" . The Bethe 
lattice, however, also has some unphysical properties. For 
instance, the description of transport is controversial 10 
for the Bethe lattice; here we choose to define the square 
of the velocity of a "quasiparticle" in such a way that the 
optical sum rule is enforced 11 . There also is no apparent 
way to determine the volume for the Bethe lattice, which 
leads to the impossibility of a microscopic determination 
of the electric conductivity unit. 

We consider the Hamiltonian for the spinless Falicov- 
Kimball model 12 
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where c\ (ci) and fj (/j) are the conduction and localized 
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electron creation (annihilation) operators respectively for 
a spinless electron at site i, and U is the on-site Coulomb 
interaction strength. The hopping integral is scaled with 
the number of nearest neighbors Z so as to have a finite 
result in the limit 13 ' 14 of Z — > oo; we measure all energies 
in units of t* = 1. As shown earlier, the case of infinite 
Z provides an opportunity to take advantage of the lo- 
cality of the self-energy and to use dynamical mean field 
theory (DMFT) (for a review see Refs. 14, 15). Here we 
consider the Bethe lattice so t he non interacting density of 
states is a semicircle p(e) = \/~A — e 2 /2ir. The method for 
solving the Falicov-Kimball model using DMFT is given 
elsewhere 15-17 . We work with a fixed value w\ = (n/) 
for the average number of localized electrons, which is the 
so-called binary alloy picture. The other two input pa- 
rameters needed to calculate the Green's function G(u) 
and self energy £(w) are U and p e — the average num- 
ber of conduction electrons. All our calculations are per- 
formed for the case p e = l — w\. This case is of particular 
interest because it yields a correlated insulating state at 
values of U larger than a certain critical U c . 

The local Green's function G = G(uj) on the real axis 
satisfies the following cubic equation 18 ' 19 : 

G 3 - 2(w + n - U/2)G 2 + [1 + (w + n - U/2) 2 

-^-}G-{u + p-U/2 + U(w 1 -±)}=Q (2) 

which we solve numerically in order to determine the 
Green's function, and subsequently, the self-energy from 
the relation (valid on the Bethe lattice): 



£(w) = uj + p-G{uj) 
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G(u>) 



(3) 



Note that the physical root must be chosen in Eq. (2) 
to yield the retarded Green's function and a chemical 
potential p is employed to get the right average electron 
density p e . 

In order to find the transport, we start from the exact 
expression for the scattering time on the Bethe lattice 11 



-{o>) = \j dep(e)(4-e 2 )A 2 (e,u;), 



(4) 



to zero outside of the band for large \lu\, whereas t(u>) 
on the hypercubic lattice approaches a nonzero constant. 
There is, however, an essential difference, namely, it can 
be shown that in the vicinity of the band edges at zero 
temperature t(uj) obeys 



t{u) rts A t 9(uj 
+B t 8(-lu 
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where A T and B T are constants, 9(x) is the unit step- 
function and Eg is the bandgap. On the other hand, 
the interacting DOS [pint(w) = — -ImG(w)] at T = 
behaves as 19 
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where A and B are again constants. Note, that at fi- 
nite temperatures the expressions in Eqs. (6) and (7) are 
modified by a temperature-dependent shift in the chem- 
ical potential (u> — > u> + p — pr=a)- This difference in be- 
havior of t{uo) and pi n t{u) in Eqs. (6) and (7) arises from 
the extra velocity squared terms in t{uo) which modify 
the behavior at the band edges. Hence, r(w) ^ Topi n t{oj) 
which would be the simplest approximation. 

Once t(u>) is known, we can compute the transport co- 
efficients Lij , according to the Jonson-Mahan theorem 20 
(see also Ref. 7): 
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where <jq has units of conductivity e 2 /ha, and a is a 
length scale that cannot be independently determined 
on the Bethe lattice. Here, f(u>) = 1/[1 + exp(w/fcsT)] 
is the Fermi-Dirac distribution. Thermoelectric charac- 
teristics can be obtained once these transport coefficients 
are determined. Thus, 
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where A(e, uj) = — -Im- 



h, — is the spectral func- 
tion. This expression can be evaluated in terms of the 
local Green's function as 



S= - 



k B L 
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The scattering time computed on the Bethe lattice pos- 
sesses several important properties. It resembles the in- 
teracting density of states in shape and behavior as a 
function of the Coulomb interaction U and filling W\. 
Namely, it develops a well-defined gap at the metal- 
insulator transition, as opposed to the hypercubic lattice 
case, where r(w) assumes a power law behavior. In ad- 
dition, on the Bethe lattice, the relaxation time is equal 
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Here, Ode and K e are the electric and thermal (electronic 
part) conductivities respectively, S is the thermopower, 
ZT is the thermoelectric figure-of-merit, and C is the 
Lorenz number. Note that the lattice contribution to the 
thermal conductivity has been neglected here. 

At low temperature, in the correlated insulator, we can 
determine analytic expressions for the electronic trans- 
port. First, note that = /i — ht=o, the change in 
chemical potential with temperature, is a linear function 
to leading order 

A^«|ln(|) + 0(T$). (14) 

Let us denote 

a = Ar \[j h = Br {l3' (15) 
then the thermoelectric characteristics follow as 

a dc ^n e-^[T{a + b)+O{Ti)], (16) 



suggest that zero temperature values of ZT can be made 
as large as desired by choosing proper values of U and 
w\. This result is, however, undermined by the fact that 
conductivity is exponentially small for low temperatures 
[Eq. (16)], so high voltage would be needed to operate a 
device. Note that on the hypercubic lattice 8 ZT —* as 
T — > in both the metal and the insulator. 
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FIG. 1. Thermoelectric figure-of-merit as a function of U 
at T = for the Falicov- Kimball model at p e — 1 — wi, 
and wi = 0.05 (solid), 0.07 (dashed), 0.1 (dotted), 0.2 
(chain-dotted), and 0.3 (long-dashed). 



and 

ZTxs ( g ~f) 2 + 0(T*). (20) 

Therefore, at low temperatures, the thermopower S and 
Lorenz number C diverge as 1/T and 1/T 2 , respectively. 
Then, Ode and n e assume exponentially small values, and 
ZT approaches a nonzero constant (in the metallic phase 
ZT vanishes as T 2 ). 

A divergence of the thermopower does not necessarily 
lead to unphysical consequences. For instance, the Peltier 
heat between two materials satisfies Q p = —T(S2 — Si)j, 
with j the electrical current flowing from material 1 to 
material 2. Even if the Peltier coefficient — T(S , 2 — Si) 
is finite as T — ► the Peltier heat vanishes because the 
current decreases exponentially with T as T — > 0. 

Now we show our numerical results. In Fig. 1 we 
present the results for the thermoelectric figure-of-merit 
ZT at T — as a function of U for several different 
values of ui\. For small U (in a metal) ZT vanishes as 
T — > 0. As the gap opens, ZT starts to grow linearly 
with U. It is clear that the asymmetry of the interacting 
DOS, which increases aswi differs from 0.5 brings about 
larger values of ZT for a given U. The closer w\ is to 
zero or unity the larger ZT is. Our numerical results 



Figure 2 shows the thermoelectric figure-of-merit and 
the Lorenz number as functions of temperature for U = 2. 
At half-filling U = 2 is the critical value of the Coulomb 
interaction for the metal-insulator transition and the crit- 
ical U c decreases as Wi changes from 0.5. (For other 
values of U > U c the temperature dependence of all 
thermoelectric parameters behaves similarly, with only 
quantitative differences.) The thermoelectric figure-of- 
merit grows monotonically with temperature and reaches 
large values at high temperature. The low temperature 
peak found close to half-filling on the hypercubic lattice 8 
does not occur on the Bethe lattice. The largest values 
of C are achieved at low temperatures in the insulating 
phase, while in the metallic phase, the Lorenz number 
has a maximum at finite temperatures. 
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FIG. 2. (a) Thermoelectric figure-of-merit and (b) 
Lorenz number for the Falicov-Kimball model at U = 
p e = 1 — wi, and mi = 0.5 (solid), 0.4 (dashed), 
(chain-dotted), 0.2 (dotted), and 0.1 (long-dashed). 
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FIG. 3. (a) Electric conductivity, (b) electronic contribu- 
tion to the thermal conductivity, and (c) thermopower for the 
Falicov-Kimball model at U = 2, p e = 1 — ioi, and wi = 0.5 
(solid), 0.4 (dashed), 0.3 (chain-dotted), 0.2 (dotted), and 0.1 
(long-dashed). 

The temperature dependence of the electric conductiv- 
ity <?dc, the electronic part of the thermal conductivity 
K e , and the thermopower S for U = 2 and various w\ are 
shown in Fig. 3. Both Ode and K e are exponentially small 
at low temperatures in agreement with the analytical ex- 
pressions in Eqs. (16) and (18), reach a maximum, and 
then decrease at higher temperatures. Larger maximum 
values of achieved closer to half-filling, contrary 

to the behavior of <Jdc, where larger maximum values are 
reached away from w\ = 0.5. The thermopower increases 
away from half-filling as well, and exactly vanishes for 
w\ = 0.5 (due to particle-hole symmetry). It has a pole at 
T = 0, according to Eq. (17). Also, the thermopower is 
positive, indicating hole-like transport for w\ < 0.5, and 
it satisfies the relation S(wi,U,T) = —5(1 — w\,U,T) 



from particle-hole symmetry. 

A comment is in order as to the influence of the ne- 
glected lattice component of the thermal conductivity. 
At low temperatures (T well below the Debye temper- 
ature Qd) the lattice contribution to the thermal con- 
ductivity kl decreases 21 as T 3 . At the same time, the 
electronic contribution K e is exponentially small, accord- 
ing to Eq. (18). As a consequence, the lattice component 
is likely to dominate the thermal conductivity and, there- 
fore, will reduce the value of ZT. On the other hand, at 
high temperatures (T > 9d) the rate of decline of the 
lattice component is 21 kl ~ 1/T X , where x is between 1 
and 2. The derivation of the precise behavior is rather 
complex, related to the competition between cubic and 
quartic anharmonic scattering processes 22 . The impor- 
tant point, however, is that at the same time n e assumes 
large values and should become dominant, allowing the 
high values of ZT to be experimentally achievable. 

In summary, we have analyzed the thermal transport 
properties of strongly correlated materials within the 
Falicov-Kimball model on the Bethe lattice. Significant 
values of ZT can be obtained for both low and high tem- 
perature regimes, contrary to the hypercubic lattice case, 
where ZT declines at high temperatures. On the other 
hand, K e grows linearly at high temperature on the hy- 
percubic lattice, while on the Bethe lattice it decreases. 
The behavior of the interacting density of states Pi n t{w) 
as well as that of the scattering time t(u>) suggest that the 
results obtained on the Bethe lattice are closer to what 
happens in real materials. However, while high ZT val- 
ues are necessary for practical applications, they are not 
sufficient to guarantee a viable device. In order to take 
advantage of these properties, one needs to be able to 
achieve appreciable values of the current densities. Since 
the conductivity approaches zero for T — > 0, the feasi- 
bility of thermoelectric devices based on strongly corre- 
lated materials described by this Falicov-Kimball model 
is questionable at low temperatures. At high tempera- 
tures, judging by the high values of ZT and lower lat- 
tice contribution to the thermal conductivity, it may be 
possible to employ these materials for power generation 
applications. 
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